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Abstract 

The electric dipole moments (EDMs) of charged leptons are significantly suppressed in stan- 
dard model. It has been found previously that they are even more severely suppressed in seesaw 
type models by powers of tiny neutrino masses as far as a leptonic CP source is concerned. We 
investigate whether a Majorana-type Yukawa coupling between charged leptons and a doubly 
charged scalar can contribute significantly to their EDMs. An observable EDM would then 
help unravel the Majorana nature of neutrinos by a lepton number conserving quantity. We 
find that the EDMs are indeed parametrically large, of the form da °^ ema{m^ — my)/m^ up 
to logarithms, where and m are respectively the masses of charged leptons and the scalar. 
And they satisfy a sum rule to good precision, de/mg + d^/m^ + d-c/ttiT: = 0. With the most 
stringent constraints from lepton flavor violating transitions taken into account, their values are 
still much larger than the mentioned previous results. Unfortunately, even in the most optimistic 
case the electron EDM is about three orders of magnitude below the foreseeable experimental 
sensitivity. 
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1 Introduction 



The CP violation through a Dirac phase in the CKM matrix of weak interactions has been well 
tested in the flavored systems of hadrons. It is generally believed however that this cannot be 
the unique or even the dominant source of CP violation because of the observed large baryon 
number asymmetry in our universe (BAU). One of attractive solutions to BAU is offered by the 
mechanism of leptogenesis [[H in which the lepton number asymmetry is first generated through 
CP violation in the lepton sector and then converts partly into BAU via sphaleron effects ^ . 

Our current information on the leptonic mixing matrix comes dominantly from experiments 
of neutrino oscillations [3]. The matrix is CKM-like involving a single Dirac phase if neutrinos 
are Dirac particles, but can contain additional two CP phases if neutrinos are of Majorana nature. 
The oscillations are blind to the latter Majorana phases while their sensitivity to the Dirac phase 
is seriously diminished by a very small, if not vanishing, mixing angle out of three. This leaves 
CP violation in the lepton sector largely untested so far except perhaps for the experiment of 
neutrinoless double beta (Ov2/3) decays which can be sensitive to CP phases but whose status 
is under debate. 

Nevertheless, there is another physical observable, the electric dipole moment (EDM), that 
can provide an independent probe to CP violation. The current experimental limits on the EDMs 
of the mercury atom flU, neutron [5j, electron [6] and muon |7l are already very impressive, 
and further improvements are expected to take place in the near future [8|. These EDMs can in 
principle be induced by the Dirac phase in the CKM matrix of standard model (SM). However, 
it was known long ago that the electric and chromoelectric dipole moments of quarks vanish to 
the two-loop order |l9l. This was interpreted as a joint result of two features in SM ifTOll . namely 
the unitarity of the CKM matrix and the purely left-handed chirality of the charged current, and 
relaxation of any of them would yield quark EDMs at two loops ifTTI . The lepton EDMs then 
become extremely small in SM as they are first induced at four loops [fT2|. This makes them a 
potentially ideal place to search for CP violation in the lepton sector. 

If neutrinos are Dirac particles, the EDMs of charged leptons will be hopelessly tiny. The 
case of quarks in SM repeats in the lepton sector in an even worse manner since neutrinos 
can be considered degenerate to very good precision at the weak scale, in which case there 
is effectively no CP violation in the lepton sector. But the situation could be different when 
neutrinos are Majorana particles because of peculiarities with Majorana CP phases lfT3l . Since 
Majorana phases dominate over the Dirac phase in this case, an observable lepton's EDM would 
not only discover CP violation in the lepton sector but also expose the Majorana nature of 
neutrinos by a lepton number conserving quantity in sharp contrast to Ov2/3 decays. Indeed, 
as pointed out in ffT4ll . there is a topologically new type of two-loop Feynman diagrams when 
neutrinos are Majorana particles that can contribute to the charged lepton EDMs. But it was 
found subsequently that this type of contribution is always severely suppressed by neutrino 
masses from virtual loops whether one works in the standard type I seesaw model [[TSl or one 
augmented with an additional Higgs doublet [ [T6| . or in type II seesaw [17J. The obtained 
numbers are actually even smaller than the four-loop result due to the CKM phase, and thus 
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would not be observable in any foreseeable experiments. 

When neutrinos are Majorana particles, the lepton number may be violated either by a bare 
Majorana-type mass of heavy neutrinos that are singlets of SM, or by some other fields that are 
active in SM and couple in particular to leptons. The physics at low energies is much richer 
in the latter case. And the simplest choice would be to add a scalar triplet as in the type II 
seesaw model lITSl . We are thus motivated to consider the most general Majorana-type Yukawa 
couplings of charged leptons to a doubly charged scalar. These couplings could arise as, but 
are not restricted to, part of interactions as in type II seesaw or a larger extension of SM. The 
CP violation encoded in the couplings can induce EDMs to charged leptons, and we find that 
this contribution is indeed parametrically large. Besides the product of Yukawa couplings and 
a lepton mass factor for chirality flip, the EDM is suppressed by charged lepton masses squared 
over four powers of the scalar mass and is partly enhanced by a logarithm. In particular, it 
incurs no suppression by neutrino masses since no neutrinos appear in virtual loops. This is the 
largest term in lepton EDMs, to our knowledge, coming from a flavor dependent CP source. 

The paper is organized as follows. In the next section we describe the Majorana-type 
Yukawa couplings between the charged leptons and doubly charged scalars, and count the num- 
ber of independent physical parameters. The two-loop diagrams for EDMs are then evaluated 
analytically in section 3, and a sum rule is found. Using the most stringent constraints from 
lepton flavor violating decays we estimate in section 4 the largest allowed values for the EDMs. 
We summarize and conclude in the last section. 

2 Majorana-type Yukawa couplings 

The relevant interactions for our study are the Majorana-type Yukawa couplings 

^Yuk = fb'f^Pd^ ++ + Ib^^PRp^ — (1) 

and the standard QED 

^QEJ, = -eAHY^l + lieA^{^-d^^++-^++d^^-) (2) 

Here £, t, and are respectively the charged lepton, doubly charged scalar and electromag- 
netic fields, and e is the electromagnetic coupling. We use Greek letters to denote the three 
charged leptons. = i^y^ is the matrix employed in charge conjugation and Pr^l = ^(l ± 75) 
are chiral projectors. 

The Yukawa coupling matrix b is symmetric in lepton flavors due to antisymmetry in fermion 
fields but is otherwise arbitrary. With n flavors of leptons, b has generally n-\-^n{n — \) moduli 
and n + ^n{n — 1) phases. All moduli are physical parameters. However, not all of the phases 
are physical. For instance, the phases in the diagonal entries baa can all be removed by re- 
defining complex fields £«• After this, there are no more degrees of freedom to rephase fields 
without reintroducing phases into baa- There are thus only ^n{n — 1) physical phases. They 
signal T and CP violation as we analyze below. 
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When the matrix b is real, we can prescribe the T and CP transformations as T^^^T^^ = 
+^'^~^, {CP)^^^{CP)^^ = so that both are preserved by the above interactions. If 

b is purely imaginary, we can prescribe in the opposite manner to preserve both T and CP, 
T^++T-^ = -<^++, (CP)<^++(CP)"i = +<^"". The latter case can of course be reduced to 
the former by rephasing the ^ field by a factor of i. Therefore, T and CP symmetries are 
violated only when the matrix b is genuinely complex, neither real nor purely imaginary. 

The above results are general and do not rely on any model. It is also possible to preserve 
lepton number by assigning two units to ^ . Our later calculation of EDM applies to the 
general case. But since the doubly charged scalars appear naturally in type II seesaw, it is 
interesting to consider this particular case separately: 

b=:^V*myV^ (3) 

2V3 

where V is the lepton mixing matrix and my the diagonal neutrino mass matrix with real, semi- 
positive eigenvalues m,. The vacuum expectation value of the scalar triplet, V3, is induced from 
that of the scalar doublet through a soft lepton number violating term. It is possible and common 
practice in type II seesaw to arrange order one Yukawa couplings b by assuming V3 to be the 
same order of magnitude as my. The moduli in b correspond to n neutrino masses (over |v3|) 
plus jn{n — 1) mixing angles in V, while the physical phases are equivalent to ^(n — 1) (n — 2) 
Dirac phases and (n — 1) Majorana phases in V. 



3 Evaluation of electric dipole moments 

Now we calculate the EDM da induced for the charged lepton ia due to interactions in eqs. ([H 
O. The effective EDM interaction is defined as 

^EDM = -^daiaY5<y^viaF^^ (4) 

There is no contribution at one loop level since the matrix element b^p always appears in a 
self -conjugate form, {b^p p, so that no phases can survive. The other way to see this is to notice 
that, when computing da for a specific a, one can choose suitable phases for the ip fields so 
that all of bap are real. Thus more factors of b have to be involved to induce an EDM, and the 
first contribution occurs at two loop level. 

The two-loop Feynman diagrams contributing to da are depicted in Fig. 1 . The incoming 
and outgoing momenta of the are respectively p±jq with q being the outgoing momentum 
of the photon attached at the vertex indicated by (g). The arrows in the graphs denote the flow of 
negative charges, and the summation over the virtual charged leptons ip, iy and i§ is implied. 
We find that because of the chirality structure in J^Yuk the chirality flip required by the EDM 
has to be done by the external lepton mass, m^. Upon extracting out this mass factor we ignore 
further dependence on it. This is a good approximation for practical purposes with incurred 
relative errors of order 0{ra), where ra = in^/m^ and m is the mass of The dependence 
on other charged lepton masses enters in a quadratic form, i.e., via rp y §. 
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The lepton flavor dependence in the relevant term of each graph can thus be described as 

f^Ya^dp^Srbjiaf{rii,ry;r§), (5) 

where / is a function of the indicated mass ratios. / is generally a sum of terms that are 
respectively symmetric and antisymmetric in /3 and 7. The symmetric term cannot contribute 
to EDM since we are effectively summing the self-conjugated b factors, ^ya^5/3^5y^j8a +c.c., 
which do not vanish even for a real or purely imaginary b. The antisymmetric combination on 
the other hand survives only when b is genuinely complex with CP phases involved: 

^3 [byablpbgybpcc] [f{rp , ry\ r§) - f{ry, rp;r§)] (6) 

Since rg <C 1, the leading term, if not vanishing, is obtained by setting r§ = 0. The b factors 
then degenerate into the form 

Wyabpaib'^b)p,] (7) 



/ / 



(a) (b) 



(e) 

Figure 1. Diagrams contributing to EDM of £«• 

It is interesting that the above form does not vanish even in the case of two flavors where 
only a single Majorana phase can appear. To see the point, it suffices to consider the easier case 
of type II seesaw in eq. ([3]) with 
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where c, s are the cosine and sine of the mixing angle, and u is the CP phase. Then, we find for 
instance 

{2\v3\)^2iZ[bl,bee{b''b)e^] = c^s^{mj-ml)mim2{u^-u*^) (9) 

which does not vanish in general. This is a feature pertaining to the Majorana-type couplings 
of charged leptons in eq. ([T]) or the Majorana nature of neutrinos in type II seesaw. 

We will see in the next section that the combination b^^bgy is no less constrained than 
{b^b)py. It is thus a good approximation to keep the leading term at = while ignoring small 
corrections that are at most of order rg Inrg. The final answer for da thus looks like 



da = C-^Z [b;abpaib^b)py\ [f{rp , ry, 0) - /(r^, rp'O)] (10) 

where C is a loop factor. This result entails an interesting sum rule 

dp du dt 

— + -^ + — = (11) 
rUe nxx 

which is exact up to small relative corrections of Oija^^)- And up to logarithmic enhancements, 
we have approximately, 

da-C \l^Z^J,^^[h^h)p^\ (12) 

We are now ready to present the results for the graphs. Graph (a) is symmetric in /3 and 7, 
and does not contribute to EDM. Graphs (b) and (c) each contain symmetric and antisymmetric 
terms, while the sum of (d) and (e) is antisymmetric. The contribution to EDM is 

25 ^ 

^« = 24(4^y^2 ^ ^abpa{b'b)py\J{r^,ry) (13) 
where again summation over /3 , 7 is implied and 7 is a sum over four graphs: 

J{rp , vy) = y (rp , ry) + J^^'^ [rp , ry) + 7^^+^) {rp , r,) (14) 
Each of these four graphs has ultraviolet sub-divergences. In 4 — 2£ dimensions, they are 

7(&)div ^ _27(0div ^ _2jid+e)div ^ F^rp,ry)r{£), (15) 

where the arguments in 7 are suppressed and 

2[b^ + c^-bc- bc{b + c)+ b^c^] 



F{b,c) 



{b-lY{c-lf{b-c) 
b^[-3c + b{l+b + c)]\nb c^[-3b + c{l + b + c)]\nc 



{b-l)\b-cY {c-l)\b-cf 
The divergences are canceled on summation as they must be. 
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The analytic result for the finite part is much more lengthy. In addition to the displayed 
function F , each graph involves one or two other twofold parameter integrals that can be worked 
out in terms of the fractions, logarithms Inr^g and Inry, and the dilogarithms Li2(l — rp) and 
Li2(l — Ty). We will not record these exact results but the sum of all graphs that has been 
expanded to the leading order in r^ y. 

rl+r^y-rpry rj{rp -3ry)\nrp - r^^{ry-3rp)\nry 

where the dots stand for higher order terms in rp y. Since the charged lepton masses are hierar- 
chical, further expansion is possible; for 1 ^ ^ ry, we have 

Jirp.ry) = rp - 2ry+^{rp - 3ry)\nrp + ■ ■ ■ (17) 

We have tested that the leading terms shown in eq. (fT6l) recover the first three digits of the exact 
results at m = 200 GeV and are good enough for our later numerical analysis. 



4 Numerical analysis 

Our result for the charged lepton EDMs shown in eqs. (I13I16I) is suppressed by charged lepton 
masses squared over four powers of the scalar mass, and has a mild logarithmic enhancement 
factor. This is a parametrically large contribution. For instance, at our reference point m = 
200 GeV, we have 7(rg,r^) ^ 1.83 x 10"^ 7(re,r^) ^ 7(r^,r^) ^ 2.94 x 10"^ and 

~ 4 X 10"^^ X [b factors] e cm (18) 

which would be within the reach in the next generation of experiment at the sensitivity of order 
10-31 e cm m. 

However, the same Yukawa couplings in eq. ([B induce other effects as well, and a realistic 
estimate of EDM should take into account the constraints from those effects. In this section, we 
present our numerical results in two approaches. The main constraints considered are from lep- 
ton flavor violating (LEV) decays of charged leptons. Also mentioned are anomalous magnetic 
moments and Ov2j8 decays. We start with a model independent analysis in the next subsection 
and then specialize to the case of type II seesaw. The constraints in the second approach are 
more stringent because of less free parameters involved. 



4.1 Approach 1: model independent result 

The Yukawa couplings in eq. ([T]) mediate radiative LEV decays at one loop level and purely 
leptonic decays at tree level. The branching ratio for the radiative decay is 



33 



a 



Gfirr 



(19) 
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modes 


/i ^ ey 


X ey 




/i ^ 3e 


T 3e 


Br 


1.2 10-^^ [201 


1.1 10"' [21| 


4.5 lO-'^ li22J 


1.0 10-^^ |23| 


4.3 10-« [24| 


bounds 


1.2 10"^ 


2.9 10"^ 


1.9 10^ 


2.0 10"^ 


1.0 10^ 


modes 


T ^ 3/i 




T ^ \lle 


T ^ ee/i 


T ^ /i/ie 


Br 


5.3 10-*^ Hal 


5.6 10-^ [24J 


5.8 10-*^ L24J 


8.0 10"^ [24J 


3.7 10"*^ [24| 


bounds 


1.1 10-^ 


1.1 10"^ 


1.2 10--^ 


9.7 10-^ 


6.6 10^ 



Table 1 : Experimental upper bounds on branching ratios of decays in eqs. (fT9l |20|) set upper 
bounds on |(Z7'''Z7)c(j(3|/(Girm^) and \h^f^h^y\l{GFn9-^ respectively. 



with 5^ = 1 and B^: ^ 17%. 5^ is a model parameter which equals (8/9) for the contribution 
of (^^^ alone (in this subsection) and equals 1 when both b,'^'^ and are included as in type 
II seesaw model (in the next). The branching ratio for the purely leptonic decay is 



Brf 



1 

22 



(20) 



which is only induced by exchange. The factor (2 — 5py) distinguishes between identical 
and nonidentical particles in the final state. Using the experimental bounds on the branching 
ratios we can constrain \{b^b)ccii\/iGfnP-) and \bsabjiy\/{Gfm^) respectively. These numbers 
are shown in table 1, and will be employed to set conservative upper bounds on EDMs. 

Each da has three terms proportional to J{re,r^), J{re,r^), and 7(r^,rT:) respectively, for 
instance, 



96;r4 J, 
G^P erue 



+3 
+3 
+3 



h* h 




irP-Gf 


rrP-Gf 




{b^b)e,' 


nP-Gp 


nP-Gf 




[b^b)^. 


nP-Gf 


nPGf 



(21) 



The first term is much smaller because of a smaller J factor and more severely suppressed 
moduli of the products of b factors, and can safely be ignored. In the optimistic case where the 
products of b factors in the last two terms are purely imaginary and add constructively, we get 
at m = 200 GeV, 



\dJ < 8.1 X 10"^^ ecm 



(22) 



Since our bounds in table 1 are given independently of itP while rn^J depends only logarithmi- 
cally on rrp, the above bound is stable against mild variations of m. Similarly, we obtain 



\d^\<lAx lO^^^ecm 



(23) 
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The expression for contains several combinations of b factors that cannot be constrained in 
LFV decays, so that a direct bound is not possible. But we can utilize the sum rule (fTTI) to set a 
bound 

<5.2x 10"^^ ecm (24) 

The limit on \de\, though larger than the four-loop SM result lfT2ll and the bounds reached via 
other mechanisms |[T4l[T5l[T6l[T7l, is still about three orders of magnitude below the precision 
reachable in the near future lfT9l . 



4.2 Approach 2: a case study in type II seesaw 

The discussion in the previous subsection is model independent. When the Yukawa interaction 
in eq. © is part of a complete structure in a model, more stringent constraints on EDMs can be 
obtained. This is the case particularly in the type II seesaw model where the Yukawa couplings 
are related via eq. ([3]) to the neutrino masses and mixing matrix which have been determined to 
certain extent. In this subsection we will not attempt a global fitting but demonstrate the point 
by a case study in this model. 

The mixing pattern determined by oscillation data is close to the tribimaximal texture ll25l . 
We will work in this simplified scenario. There is then no Dirac phase but there can be two 
Majorana phases u\ 2'- 



V 



( 



\ 



'2„ 
3MI 



73"2 
73"2 



M2 





1 

72 



^/2 / 



Then, the matrix h^^h is real and symmetric. 



(25) 



4|v3|VZ7 = -(mf + m^ + mf)l3 



+- 




-2Ai2 

-Ai3 
-Ai3 - 2A23 




2 



(26) 



with A/j 

The off-diagonal moduli \ {h^h)af> \ depend explicitly on A,y, which have been determined 
e.g. in [|26l to be, A21 = 7.6 x lO'^ eV^, IA31I = 2.4 x lO"^ gY2 jj^g y^oxxxidi on Br(/i ej) 
then implies (using B^ — \) 



|v3|WGf>5.75x 10^2 



(27) 



Since Br(T ey) also depends on A12, its less stringent bound is useless. Instead, its relation 
to Br(/i — i> ey) in type II seesaw and the experimental bound on the latter mean 

12 



Br(T ey) = BrBr{{i ey) < 2.0 x 10" 



(28) 
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which is much below the current bound. We also notice that the bound on Br(T jiy), though 
more than three orders of magnitude larger than Br(/i — > ey), gives a constraint that is only 
slightly weaker than in eq. (ITtI) . because of an enhancement factor IA31 1/A21. A similar relation 
also holds between Br(T — > 3e) and Br(jU 3e) ; the more stringent bound on the latter implies 

Br(T 3e) = B^Bx{n ^ 3e) < 1.7 x 10"^^ (29) 

which is much smaller than its current bound and thus more difficult to observe than other decay 
modes of T. 

To proceed further with leptonic decays of x and EDM, we set all neutrino masses in b*g^bpy 
(but not in {b'^b)pcc of course) to be equal to their average value my This simplification holds 
true barring very delicate cancellation among neutrino mass differences and Majorana phases. 
Then, the most stringent constraints on /i ^ 3e, T ^ ee/i, fLjle (together with the less stringent 
one on T ^ e2/i) are proportional to |mj — Wjl- We may reach the most optimistic values of 
EDMs by assuming = = e"^ to avoid these bounds. The remaining ones on T ^ 3/i, file 
yield comparable constraints: 



-2 

my 


sin0 




8| 


\V3\ 





3 


|sin(<^)/2)| 


' 4| 


V3 


\^rrP-Gf 



xlO-3, ""7"^ I < 1.2x10-3 (30) 



The electron EDM being proportional to Ziuf-u^) vanishes, while the other two simplify to 

di m^Ai3sin^ 



Jir^i^r,), (31) 



barring cancellation of 0(A2i/|A3i|) ~ 3% or 0([A2i7(r„ r^)]/[|A3i|7(r^, r^)]) ~ 2 x 10-^. 
The bounds in eqs. (l27l l30l) then give at m = 200 GeV 



\d^\ < 2.0 X 10-33 e cm, \dr\ < 3.4 x lO^^^ g cm (32) 

As expected, this result is better than the model-independent one in the previous subsection. 

In this special scenario, the effective neutrino mass measured in Ov2/3 decays is ifij^j^ — 
(2mi +m2)/3 ~ my The contributions to anomalous magnetic dipole moments depend on the 
diagonal elements of b^'b and are given by 

- - ^ "^'"'^^ < 1.1 X 10-14 



(4;r)2 4|v3|2m2GF 
3 mlm}.Gf 

a = ,9 9^ < 4.7 X 10-1° (33) 

^ (4;r)2 4|v3|2m2GF 

using eq.(|27]) and rhy ~ 0.21 eV from a recent update of cosmological bounds on the sum of 
neutrino masses ETl . Both ag and are below the potential gap between measurements and 
SM expectations Il28ll29ll. 
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5 Conclusion 



CP violation in the lepton sector has remained an experimentally unexplored issue. The charged 
lepton EDMs offer a potential arena to detect it. This is especially encouraged by the experi- 
mental precision in EDM measurements that has been reached and will possibly be accessible. 
However, it has been found previously that it is hard to obtain a not too tiny EDM for charged 
leptons from a flavor CP source. This may be blamed on the very light, thus almost degenerate 
at the electroweak scale, neutrinos. Together with a small mixing angle out of three, this makes 
a Dirac phase effectively unobservable; and it suppresses the effects of Majorana phases on 
EDMs by several factors of neutrino masses. We have thus been motivated to consider a CP 
source that arises from Majorana-type Yukawa couplings of charged leptons. Such couplings 
may appear naturally in SM with an extended scalar sector, such as type II seesaw model, but 
we have presented our analytic results in a general setting. We found that the EDMs so obtained 
are parametrically large. They are only suppressed by charged lepton masses squared over four 
powers of heavy scalar masses for order one Yukawa couplings that may be naturally arranged, 
for instance, in type II seesaw model by assigning a tiny vacuum expectation value for the scalar 
triplet. 

Nevertheless, the fate with a flavor CP source seems insurmountable. While a large enough 
EDM, though flavor diagonal, demands reasonably large flavor changing couplings, this may 
not be allowed by strictly bounded LEV transitions. With these bounds taken into account, we 
found that the electron EDM is at least three orders of magnitude below the precision achievable 
in the near future, although it is still much larger than the contributions considered previously. 

Acknowledgement This work is supported in part by the grants NCET-06-0211 and NSFC- 
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